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Abstract 

Hollerbach and Riidiger have reported a new type of magnetorotational instability (MRI) in 
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\^ • magnetized Taylor-Couette flow in the presence of combined axial and azimuthal magnetic fields. 
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, The salient advantage of this "helical" MRI (HMRI) is that marginal instability occurs at arbitrarily 
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low magnetic Reynolds and Lundquist numbers, suggesting that HMRI might be easier to realize 
than standard MRI (axial field only). We confirm their results, calculate HMRI growth rates, 
and show that in the resistive limit, HMRI is a weakly destabilized inertial oscillation propagating 
in a unique direction along the axis. But we report other features of HMRI that make it less 
^ - attractive for experiments and for resistive astrophysical disks. Growth rates are small and require 

^ . large axial currents. More fundamentally, instability of highly resistive flow is peculiar to inflnitely 

long or periodic cylinders: flnite cylinders with insulating endcaps are shown to be stable in this 
limit. Also, keplerian rotation proflles are stable in the resistive limit regardless of axial boundary 
conditions. Nevertheless, the addition of toroidal fleld lowers thresholds for instability even in flnite 
cylinders. 
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I. INTRODUCTION 



The magnetorotational instability (MRI) is probably the main source of turbulence and 
accretion in sufficiently ionized astrophysical disks (l| . MRI was ffist discovered theoretically 
RO; then later supported numerically 0, ^,0], but has never been directly observed in 
astronomy. No laboratory study of MRI has been completed except for that of Sisan et al. , 
whose experiment proceeded from a background state that was not in MHD equilibrium. We 
and others therefore have proposed experimental demonstrations of MRI 00, Q. The 
experimental geometry planned by most groups is a magnetized Taylor- Couette flow: an 
incompressible liquid metal confined between concentric rotating cylinders, with an imposed 
background magnetic field sustained by currents external to the fluid. 

The challenge for experimentation, however, is that liquid-metal flows are very far from 
ideal on laboratory scales. While the fluid Reynolds number Re = Qiri{r2 — ri)/v can be 
large, the corresponding magnetic Reynolds number i?em = f2iri(r2 — ri)/// is modest or 
small, because the magnetic Prandtl number Pr^a = vjri ^ 10^^ — 10"^ in liquid metals; 
here v < 10^^ cm^ s~^ is the kinematic viscosity and r/ is the magnetic diffusivity. Standard 
MRI modes will not grow unless both the rotation period and the Alfven crossing time are 
shorter than the timescale for magnetic diffusion. This requires both i?em ^ 1 and S* > 1, 
where S' = V^(r2 — ri)/// is the Lundquist number, and Va = B / ^/ioP is the Alfven speed. 
Therefore, Re > 10^ and fields of several kilogauss must typically be achieved. 

Recently, Hollerbach and collaborators have discovered that MRI-like modes may grow 
at much reduced -Rcm and S in the presence of a helical background field, a current-free 
combination of axial and toroidal field ^ Q| . 

= i?f (e. + p'j-ee) (1) 

in cylindrical coordinates {r,6,z), where Bi^^ and /3 are constants. (When it will not cause 
ambiguity, we will omit the superscript (0) from B and Bz hereafter.) Henceforth, "standard 
MRI" (SMRI) will refer to cases where the P = 0, and "helical MRI" (HMRI) to modes 
that require (3^0. In centrifugally stable flows — meaning that d{r'^Q)^ / dr > 0, where 
Q = Vg^^ /r is the background angular velocity — SMRI exists only when Re^ and S exceed 



thresholds of order unity [2l,lU|- Remarkably, however, HMRI may persist in such flows even 
as both parameters tend to zero, though not independently: more precisely, the thresholds 
are <C 1 and would vanish if the fluid were inviscid {u = 0). In a fixed geometry and fiow 



profile, the resistive limit may be approached theoretically by increasing r] with all other 
parameters held constant. The growth rate of inviscid HMRI is then oc so that the 
hydrodynamic case is a ppr oached continuously. The special case of toroidal-only magnetic 
field (/3 = oo) is stable IJ]. 

A novel feature of the background state for HMRI is that there is a uniform axial flux 
of angular momentum carried by the field, rT^^^^ = —rBgBz/ fio and an associated axial 
Poynting flux Q times this. In an infinite or periodic cylinder, the question of the sources and 
sinks of these axial fluxes need not arise, but in an experimental device, a torque is exerted 
by the axial field on the radial sections of the coil that complete the circuit containing the 
axial current. Related to this perhaps, the dispersion relation for linear modes is sensitive 
to the sign of the axial wavenumber (kz), and the instabilities of axially infinite or periodic 
cylinders are travelling rather than standing waves, as noted by Knobloch jisl . Ilq . This 
begs the question what should happen to the modes in finite cylinders, a question that has 
motivated much of our analysis. 

Even the analysis for periodic cylinders implies two practical difficulties for an HMRI 
experiment. First, as will be seen, the typical growth rates tend to be smaller than those 
of SMRI except in regimes where SMRI would also be unstable. Secondly, the axial current 
needed for the required toroidal fields tend to be quite large: /[kA] = 5Bgr[kG-cm]. 

In Section |H] we analyze the linear stability of HMRI using complementary approxima- 
tions, some for infinite/periodic cylinders and others for finite ones. The results are compared 
with one another and with fully nonlinear axisymmetric simulations. Our conclusions are 
summarized in Section UTTl 



II. LINEAR THEORY 



All magnetic fields are expressed as Alfven speeds, in other words, units such that /iq = 
1/p are used. Upper-case letters are used for the background magnetic field ((T)) and velocity 
V = rQ{r)e0, and lower-case {b,v) for perturbations. Frequently occurring derivatives are 
abbreviated by dl = dj. + r^^, D = drdl + d'^. Incompressibility allows the use of stream 
functions for the poloidal components: Vr = dzcf), Vz = ~dl(f), = dzip, bz = —dlip; note 
that these definitions differ by factors of r from the usual ones. The linearized inviscid MHD 
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equations then become, since Bz and rBe are constant, 

{d^-r]D)^ = Bzdz<j), (2) 
{dt - r,D)he = dz (^^0 + B-^ve + rO^ , (3) 

dtD<P - 2QdzVe = BADi; - ^9,6,, (4) 

r 

dtVe + r~\r^Qydz(t> = BAhe. (5) 
The underhned terms above are neffhffible in the resistive hmit, where h scales oc ii^^ com- 

n n 

pared to v. Neglecting these terms has been shown to suppress SMRI [111, im, but not 
HMRI as will be seen. 

Taking another time derivative of and eliminating dtVg via © yields 

[dlD + K'dl) ct> = BAdtDi; + 2 (nB^d^ - ^d^d^j be , (6) 

in which = r^^d{r'^Q)'^ /dr"^ is the square of the epicyclic frequency. As t] oo, (0) 
reduces to 

{drdl + dl)d^<P +n'{r)dl<P = 0. (7) 
A. WKB for infinite or periodic cylinders 

If we take the gap to be narrow, d = r2 — ri r, then it is reasonable to treat r, Bz, ^l, 
rQ' = 2RoQ, and r^^(r^i7)' = 2(1 + Ro)Q = /2VL as constants, and to look for perturba- 
tions oc exp{ikrr + ikzZ — iut) . The Rossby number Ro = ^dln^l/dlnr has been introduced. 
In this case one expects to have WKB solutions with D replaced by —{k'^ + k'^) = K^, where 
the total wavenumber K = 0{d~^). 

When applied to (|7j) {i.e. for t] oo) these prescriptions yield the dispersion relation for 
hydrodynamic inertial oscillations (hereafter 10), 

-io = -^^ where = ^^^ir^nf = 4(1 + Ro)n^ (8) 

10 exist only in the Rayleigh- stable regime > 0, Ro > —1, and their frequencies lie 
between and n. 

HMRI occurs at finite r] when Be/r = ug is comparable to kzBz = ujz- Define tu^ = rjK"^ 
and = kz/\K\ G [—1,1]. The dispersion relation corresponding to the system (0)-© is 



then 



r o o ^. ~i 

+ [uj'l^'^K^ - 4iuJr,i^guj^fi'^n{2 + Ro) + cu^ + 4/i2cc'^fi2i?o] , (9) 

where the complex growth rate s = —iu has been used so that the coefficients are all real 
except for those linear in uq. It is instructive to consider the limit in which o;^ is much larger 
than all of the other frequencies, including uj: 

+ ujjo + 2uj-^ [s^^ + (2/i2cj2 + ^2 ^ ^2^^^ _ 2iujsuj-yn{2 + Ro)] ^ 0{uj-^). (10) 
The replacement /x^k^ — > cUjq emphasizes that u) ~ icujo in this limit. The roots are 

u ^ Tcuio + [±2tU0CU^tufoV^fi(2 + i?o) - {2h'^ujI + cu^)] + 0(6^-^), (11) 

the bivalent signs being correlated. The other two roots of (0) represent rapidly decaying 
magnetic perturbations, s ~ —uj^^. 

We conclude that in highly resistive flow, HMRI reduces to a weakly destabilized inertial 
oscillation. In the present inviscid approximation, instability persists to arbitrarily large 
resistivity, though with reduced growth rate. Furthermore, we note from (fTT|l that instability 
[i.e. '^{uj) > 0] occurs only if the bivalent signs are chosen so that QBgBzkz/^{uj) < 0, which 
implies that the unstable mode propagates axially with the same sense as the background 
Poynting flux. [From (jSJ, the group velocity d^{uj) / dk^ and phase velocity 3?(u;) / k^ have the 
same sign.] Although we have derived this propagation rule in the resistive limit, numerical 
evidence indicates that it is true of the full dispersion relation ©, as demonstrated by 
Figure. [T] 

Instability requires the square brackets in pO|) to be positive, whence 

2{^ujef ± 1"^^^ ^^(/icue) < 0. 

Vl + Ro 

The inequality is possible if and only if the discriminant of the lefthand side, regarded as a 
quadratic equation in /icug, is positive: 

(i±^.;-8.;>o. 

l + Ro 



5 



Figure 1: Selected roots of full dispersion relation ® for t] = 2,000cm^s^^ [gallium], ri = 9cm, 
r2 = 11cm, vertical periodicity 2h = 16 cm, ili = 100 rpm, = 68.1 rpm, = 500 G, Bg = 10 kG 
at r = (ri + r2)/2. The two rapidly damped modes are omitted. 
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(a) Growth rate 7 — "^uj vs. Wave number kz 



(b) Real frequency ujr — ^to vs. Wave Number 



which translates to 



Ro < 2(1 - V2) ^ -0.8284 or Ro > 2(1 + V2) ^ 4.8284 



(12) 



Thus, at least in WKB, the keplerian value Ro = —3/4 is excluded, as of course is uniform 
rotation (Ro = 0). We say "of course" because, the background being current free, the only 
source of free energy is the shear. 



B. Numerical results for wide gaps in periodic cylinders 

We have adapted a code developed by [ill ^'^ allow for a helical field. Vertical periodicity 
is assumed, but the radial equations are solved directly by finite differences with perfectly 
conducting boundary conditions. The underhned terms in eqs. (0)-® are retained, and 
viscous terms are added although their infiuence is small at Reynolds numbers of interest. 
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13|. Table m compares the 



The code reproduces pubhshed results for marginal stability 
predictions of the WKB dispersion relation Q with those of this code (labeled "Global"). 
The agreement is reasonably good, considering the crudeness of the WKB approximation. 
No unstable modes are found for the parameters of Figure Q at Ro{ri) > —0.821: the 
keplerian value Ro = —0.75 is stable. 

Astrophysical disks correspond to very wide gaps, r2 — ri ^ h, as well as keplerian 
rotation. Given {Re^n, S) = (0.1,0.03) and r2/ri = 2.0, 2.83, 5.0, the maximum unstable 
Rossby numbers at the inner cylinder are found to be Ro{ri) = —0.928, —0.947, and —0.981, 
respectively, from our radially global linear code. We conjecture that keplerian flows — more 
precisely, flows in which > Ro > —3/4 at all radii — are stable for all gap widths. It would 
be interesting to prove this. 

We have also estimated a few growth rates with our nonlinear, compressible non-ideal 
MHD code which is a modified version of the astrophysical code ZEUS2D 3|. In this 
case, we use the wide-gap geometry of the Princeton MRI experiment lll|, except that 
the computation uses periodic vertical boundaries: ri = 7.1cm, r2 = 20.3 cm, h = 27.9 cm, 
Qi = 400 rpm, ^2 = 53.3 rpm, Bz = 500 G, Be{ri) = 1 kG; the material properties are 
again based on gallium: rj ^ 2000cm^s~^, u ^ 3 x lO^^cm^s"^. The growth rate and real 
frequency from the ZEUS2D simulations are respectively 1.06 s~^ and 3.93 s~^, compared to 
1.05 s~^ and 3.89 s^^ from the linear code. WKB yields (7,0;^) = (0.41, 3.90) s~^ — not an 
accurate result for the growth rate, but considering the width of the gap, the agreement is 
pleasing. 

The growth rates in Table IJ are of order ls~^, as compared to ~ 30 s^^ for SMRI in 
this geometry at the full rotation rate and field planned for the Princeton experiment jisl ]: 
Qi = 4, 000 rpm, ^2 = 533 rpm, Bz = 5kG, and Bg = 0. We have found this to be typical: 
the growth rates of HMRI are small, at least in regimes where SMRI is stable. 



C. Finite cylinders: a perturbative approach 

In finite nonperiodic cylinders with insulating or partially insulating endcaps, the MHD 
eigenfunctions are intrinsically two dimensional: they are not separable in r and z. (Sepa- 
rability could be achieved with perfectly conducting endcaps, but then the axial field would 
be attached to them. This would allow the boundary to exert magnetic forces on the fluid. 
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Table I: Comparison between WKB and numerical growth rates in a vertically periodic Couette flow 
with the parameters of Figure ^ except for a nonzero viscosity like that of gallium: u = 3.1x 10^^. 
The mode number n = kzh/n. 



n 


WKB 7 [s-i] 


WKB LOr [S-^] 


Global 7 [s-i] 


Global LOr [s ^] 


1 


0.1612 


0.9443 


0.0965 


1.4004 


2 


0.3911 


1.9182 


0.3465 


2.5164 


3 


0.5878 


2.7084 


0.6031 


3.2638 


4 


0.7387 


3.2646 


0.7907 


3.7094 


5 


0.8356 


3.6221 


0.8960 


3.9549 


6 


0.8805 


3.8366 


0.9339 


4.0799 


7 


0.8829 


3.9565 


0.9241 


4.1352 


8 


0.8543 


4.0166 


0.8831 


4.1512 


9 


0.8049 


4.0400 


0.8227 


4.1451 



which seems undesirable and in any case is experimentally less realistic than insulating end- 
caps.) The purely hydrodynamic problem for rj = oo is separable, however, if viscosity is 
neglected so that we may assume no-slip boundary conditions. This suggests a perturbative 
expansion of the eigenvalue problem in ri~^ — more properly, {Rem, S) {eRera,€S), with 
e a small parameter. The cylinders themselves are assumed infinitely long and perfectly 
conducting; although this is not realistic, it does not result in any attachment of the field 
to the boundaries, and it allows the magnetic field more easily to be matched onto vacuum 
solutions that decay as |z| ^ oo in the regions above and below the fluid. The underlined 
terms in equations (0)-© will be neglected because they contribute to the eigenfrequency 
only at 0{ri~'^) and higher orders. 

We begin with the zeroth-order problem, i.e. for rj = oo. As noted above, the hydrody- 
namic boundary conditions 

= on r = ri, r2 and on ^ = 0, h, (13) 
and inertial-mode equation ((7j) are separable, so we look for an eigenmode of the form 

(j)(t, r, z) = e"*'^V('^) sin kz, k = n- . (14) 
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The radial function y9(r) satisfies 

d?Lp 1 dip 
dr"^ r dr 



4abA;^ 



- 1 



V = 0, 



(15) 



assuming a Couette profile Q{r) = a + br~^ so that = 4afi, which is satisfied by the 
Bessel functions Jy{pr) & Yy{pr) if 



= 1 



4abfc^ 



2 ' 



(16) 



We may thus solve this problem exactly. However, for qualitative information, we notice 
that if we multiply (fTH|) by r it becomes 



dr \ dr J 



1 n,n 4abP 

— 4a^A;V H 

tu^ V r 



- + /c^r 
r 



This is the same form as the Sturm-Liouville problem 

|:(w^)+[Ai?(r)-Q(r)]^ = 



V9(ri) = 0, ip{r2) = 0, 



where 



= 0. 



Pfr) 



i?(r) = 4a^k^r + 



4abfc^ 



> 0, 



g(r) = A;V + - > 0, 
r 

A = 1/lu\ 



Therefore A is real and positive (j20|. Chapter X]); consequently the frequencies uj which we 
seek are all real. There are no modes which grow in time. Thus we conclude that, all inviscid 
axisymmetric modes are neutrally stable in the limit of infinite resistivity. The coefficient 
R{r) = $(r), the Rayleigh discriminant, so this result is to be expected. 
We may arrange for 0(ri) = by taking 



'■Pmnir) = Ju{pri)Yy{pr) - Y^{pri)Jy{pr). 



(17) 



Since we also require (pij^) = 0, the determinant 



A(a;,/c) = Ju{pri)Y^{pr2) - Ju{pr2)Yu{pri) 



(18) 



must vanish. The condition A = defines a discrete set of eigenfrequencies ooi^n > ^2,n > 
. . . > Umn . . . > for each k = kn- Let 0m,n be the complete eigenfunction (fT^ corresponding 
to a given kn & i-Om,n- We define an inner product [here (pmn is defined by (fT^ with ip{r) — >■ 

h r2 

{(t^m'n'Amn) = j^^j l^dr (f)m'n'(l)mn, (19) 
ri 

where the overbar denotes complex conjugation. The eigenfunctions are orthogonal in the 
sense that {(pmn, /tVm'n') = if u;^„ 7^ ^m'„'- 

To get the 0(?7^^) corrections to ujmn, we must express the magnetic perturbations ip 
and bo appearing on the righthand of © in terms of the zeroth-order eigenfunctions (pmn- 
Neglecting the time derivative in (j2)) yields 

Dij = -r]-'BA<Pmn. (20) 
To get bg from Q, we first use ^ to write vg ^ {2a/iu!mn)dz(pmn, so that 

bg ^ -2r,-'D;r' { — d-Amn + ^^^^n.n') • (21) 
V r ^ran J 

Note that we have replaced d1 with — A;^; we may similarly replace any even power of dz but 
not an odd power, which changes a sin k„,z to a multiple of cos knZ. The operator D^^ is the 
inverse of D with the boundary conditions appropriate to bg, which are different from those 
of [eq. (HI]: 

dlbg = at r = ri, r2 and 6e = at 2; = 0, /;.. (22) 
Using & (pT|l to ehminate D%Ij and fee from (jHl) results in 



UJrn.n I XT Ul, 



(23) 



On the righthand side of ()23|). the eigenmode and eigenfrequency have been evaluated 
to zeroth order in ri~^. On the lefthand side, we must consider that u ujmn + and 
4'mn + S4>, where 600 and 6(f) are of first order in rj'^. We may obtain an expression 
for 6uJ by taking the inner product of (j2Sl) with 0m,n and replacing idt Umn + 6uj on 
the lefthand side. The single term involving Sep at 0{rj^^) is (0mn, (/^^ — ^mn^)^^) y ^^'^ 
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this vanishes upon integration by parts. On the right side, it is convenient to define the 
self-adjoint operator 

H ^ 2 {-^id, + = H\ (24) 



r L0„ 



At last, then. 



i 



2r] 



2 1 2hBzk'^ 2 1 



■ (25) 



Now D and D^^ are negative-definite operators. Therefore, the only term that can make a 
positive contribution to the growth rate 3(5c<j) is the last term on the righthand side, and 
specifically the part of H involving Bgdz since ab > 0. 

To evaluate 6u! from (|25|). we need explicit expressions for D and D^^. The first is easy 
enough: it follows from ((Tj) that Dcpmn = ~(^n'^^('^)/'^mn)0"m- For -Dj^^, we construct the 
eigenfunctions of D with the boundary conditions 



Dxjnir, z) = -{q] + A;Dx,n(r, z), (26) 

TT 

Xjn(r, z) = Rjnir) sin knZ , ^ (^'^) 



Mqjri)Yi{qjr) -Yo{qjri)Ji{qjr), if ^ 0; 
where Rjn = \ (28) 

if go = 0; 

and qj satisfies Jo{qjri)Yo{qjr2) - Yo{qjri)Jo{qjr2) = (29) 

When applied to Xjn, {ij + arbitrary function /(r, z) can be expanded in 

these eigenfunctions, so that 

„ j VCjni Xjn) 

The important point is that turns a function proportional to sin k^z into another 

such. Therefore, (0mn, D^^dz(t>mn) = 0, and so the part of H involving r^^Boidz does 
not contribute to the expression ()25|) for the first-order eigenfrequency. This, however, was 
the only term that might have made for a positive growth rate. We conclude that at 0{r]~^), 
HMRI does not grow in finite cylinders with insulating endcaps. 
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The same perturbative method could have been used for periodic vertical boundary con- 
ditions; (pmn and Xjn would have involved exp{ikzz) instead of sin knZ. The term involving 
r'^Boidz in eq. (|25|) would then have contributed to the growth rate with the same sign as 
— {kz/uJmn)^BgBz. Evaluating this term, we conclude that in highly resistive periodic flows, 
(i) unstable modes propagate axially in the direction of the background Poynting flux — as 
found in WKB; and (ii) the instability occurs only if (3 > akzr/uOmn because of the second 
term in H. 

We have written MATLAB procedures to evaluate eq. ()25p. The results confirm our 
conclusions above. When periodic boundary conditions are used, the perturbative result 
matches the growth rate found from our radially global linear code to three digits in suffi- 
ciently resistive cases: e.g., 7 = 1.89 x IQ^^VLi in the Princeton geometry with i?em = 0.1, 
S = 0.043, ^2/^1 = 0.1325, /3 = 2. But when insulating endcaps are imposed, the pertur- 
bative estimate of the growth rate is always negative. 

D. Finite cylinders: two other approaches 

Here we analyze finite cylinders by approximations that do not require large resistivity: 
by a variant of WKB, and by direct axisymmetric numerical simulations. 

In the modified WKB approach, perturbations are again assumed to vary as exp{ikr + st) 
with a common complex growth rate s = —iuj and radial wavenumber kr = vr/(r2 — ri), but 
the vertical dependence is treated differently. With the t and r dependence factored out, the 
linearized equations of motion reduce to homogeneous ordinary differential equations with 
coefficients independent of z. Elementary solutions of these equations exist with exponential 
dependence on z; however, since the vertical boundaries are not translationally invariant, 
the wavenumber kz need not be real, and growing modes can be linear combinations of the 
elementary exponential solutions with the same uo but different kz- The vertical magnetic 
boundary conditions require the fields to match onto a vacuum solutions that decay expo- 
nentially as 1^1 00 in the space ri < r < r2 between the extended conducting cylinders: 

z = Q : (j) = bg = 0, dzip = \kr\ip; z = h : (f) = be = 0, dzi^ = —\kr\^p. (31) 

We search iteratively for such modes as follows. Given a trial value for s, the dispersion 
relation has six roots — in general complex — for the vertical wavenumber, which can be 
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regarded as algebraic functions of the growth rate: {^2,0(5)}, a G {1,...,6}. We seek a 
mode in the finite cyhnder of the form 

6 

q{t, r, z) = [0, ve, ^, hf = 6''+""^' ^ Y^q^ exp{tK,^z). (32) 

Each term in the sum above is the elementary solution corresponding to a particular root 
kz,a{s), with a 4-component column vector; these elementary solutions are superposed 
with constant weights {Ya}. Substitution into the boundary conditions ()31|) yields a sixth- 
order homogeneous linear system for the {Ya}. Nontrivial solutions exist only if the de- 
terminant T'(s) of this system vanishes. The equation T'(s) = is transcendental and we 
cannot solve it analytically, but a numerical nonlinear zero-finding algorithm recovers the 
roots for s. 

We have checked this procedure by replacing (jHTj) with periodic boundary conditions and 
comparing the results with direct solutions of the dispersion relation Q. Also, we find 
reasonably good agreement with growth rates determined from ZEUS2D simulations of a 
narrow-gap configuration with insulating boundaries (see below). However, for sufficiently 
large resistivity, no roots with positive 9ft(s) are found, in agreement with the perturbative 
results of Section III (^1 

For the ZEUS2D simulations, we represent the poloidal magnetic field at 2; < and z > h 
by fiux functions $±(r, z) satisfying brSr + b^e^ = r~^eg x V$ and V X 6 = 0. The latter 
implies rdr{r~^dr^) + = 0, which is solvable by separation of variables since we require 
$ = on the vertically extended conducting cylinders. The elementary solutions are 

$fc(r,z) oc re-^l^-^"l [Yi{kri)Ji{kr) - Ji{kn)Yi{kr)] , 

for an infinite discrete set of nonnegative values of k determined by ^k{f2,z) = 0. At each 
endcap, we match the vertical field protruding from the fiuid with a superposition of 
vacuum solutions of this form, and thereby obtain a boundary condition relating bz and br- 
Of course bg = at these boundaries since the current along the axis is constant. 

We have performed simulations with insulating endcaps for the parameters of Figure 
We find a complex growth rate s ~ 0.51 + 4.18^ s~^, as compared to s ~ 0.37 + 3.68z s~^ 
from the modified WKB approach (jS5)"(IS2I) above. Considering the approximate nature 
of the latter approach, the agreement is satisfactory. We have also carried out ZEUS2D 
simulations with insulating endcaps in the wide-gap experimental geometry [(ri,r2,/i) = 
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(7.1, 20.3, 28) cm]. Here we find a growth rate ~ 0.27s~^, as opposed to ~ 1.06 s~^ with 
periodic boundaries. We conclude that insulating endcaps lower the growth rate, even in 
flows of moderate (-Rcm, S). 

A limitation of our direct simulations is that since we use explicit time stepping, we cannot 
explore very large resistivities [18j. The modifled WKB approach does not suffer from any 
restriction on t], but it is not trustworthy for wide gaps. The concordance between the two 
approaches where both are applicable — namely for narrow gaps and moderate (-Rcm, S) — 
inclines us to trust results obtained from one of these approaches in regimes where the other 
is not applicable. In particular, the modifled WKB method predicts that highly resistive 
flows are completely stable in flnite cylinders, at least for narrow gaps. The perturbative 
analysis of Section 111 CI reaches the same conclusion for gaps of any width, but that analysis 
is valid at 0(?7~^) only. 

III. CONCLUSIONS 

We have analyzed the linear development of helical magnetorotational instability in a 
non-ideal magnetohydrodynamic Taylor-Couette flow, paying particular attention to the 
effects of the axial boundary conditions. A number of complementary approximations and 
numerical methods have been used. 

For inflnitely long or periodic cylinders, we conflrm that there is an axisymmetric MHD 
instability that persists to smaller magnetic Reynolds number and Lundquist number in the 
presence of both axial and toroidal background magnetic fleld than the standard MRI that 
exists for axial fleld alone. The new mode is an overstability and propagates axially in the 
direction of the background Poynting flux —rQEgB^/ fiQ. In highly resistive flows, the new 
mode is a weakly destabilized hydrodynamic inertial oscillation. Growth depends also on 
the ratio of shear to rotation, i.e. Rossby number: for all aspect ratios r2/ri that we have 
explored, and certainly for narrow gaps, the keplerian Rossby number is stable. Even for 
those proflles that permit growth, the rate tends to be rather small, except in flows that are 
sufficiently ideal to permit growth of standard MRI (axial field only). 

We have also considered finite cylinders with insulating endcaps, which are closer to 
experimental reality but which do not permit traveling modes that propagate indefinitely 
along the axis. Astrophysical disks also, of course, have limited vertical thickness. These 
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boundary conditions reduce the growth rate of the hehcal mode, and stabihze highly resistive 
flows entirely. The small growth rate may make it diflicult to detect the instability in the 
face of Ekman circulation and other experimental imperfections. In addition, the new mode 
requires toroidal flelds at least as large as the axial fleld, and therefore large axial currents, 
which introduces additional engineering challenges. 

For all of these reasons, the experimental advantage of helical MRI over standard MRI 
is open to question, as is the relevance of HMRI to astrophysical disks, although it may be 
relevant to stellar interiors and jets, where the magnetic geometry and the Rossby number 
may be more favorable. Also, HMRI may have theoretical signiflcance that goes beyond its 
direct applications. It is not understood why linearly and axisymmetrically stable rotating 



flows are often also nonlinearly and nonaxisymmetrical 
ical transition does occur at some Rossby numbers 



y unstable, especially since subcrit- 



2l\ . The fact that even a very poorly 



coupled magnetic fleld can sometimes linearly destabilize such flows hints that it might also 
affect nonlinear transition. 
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